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Abstract 

The photoproduction of beauty quarks in events with two jets and a muon has 
been measured with the ZEUS detector at HERA using an integrated luminos- 
ity of 110 pb~^. The fraction of jets containing b quarks was extracted from 
the transverse momentum distribution of the muon relative to the closest jet. 
Differential cross sections for beauty production as a function of the transverse 
momentum and pseudorapidity of the muon, of the associated jet and of x}^^^, 
the fraction of the photon's momentum participating in the hard process, are 
compared with MC models and QCD predictions made at next-to-leading order. 
The latter give a good description of the data. 
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1 Introduction 



The production of beauty quarks in ep collisions at HERA is a stringent test for perturba- 
tive Quantum Cliromodynamics (QCD) since the large 6-quark mass (mb ~ 5 GeV) pro- 
vides a hard scale that should ensure reliable predictions. When Q^, the negative squared 
four-momentum exchanged at the electron vertex, is small the reaction ep — e'hhX can 
be considered as a photoproduction process in which a quasi-real photon, emitted by the 
incoming electron, interacts with the proton. 

For 6-quark transverse momenta comparable to the 6-quark mass, next-to-leading-order 
(NLO) QCD calculations in which the h quark is generated dynamically are expected to 
provide accurate predictions for h photoproduction [1]. The corresponding leading-order 
(LO) QCD processes are the direct-photon process, in which the quasi-real photon enters 
directly in the hard interaction, 7(7 bb, and the resolved-photon process, in which the 
photon acts as a source of partons that take part in the hard interaction {gg — s> bb or 
qq — > bb). 

The beauty- production cross section has been measured in p'p" collisions at the ISR [2], 
SppS [3] and Tevatron colliders [4], in 77 interactions at LEP [5] and in fixed-target vrA^ [6] 
and pN [7] experiments. Apart from the SppS data and the fixed-target experiments, the 
results were significantly above the NLO QCD prediction. The HI measurement in ep 
interactions at HERA [8] found a cross section significantly larger than the prediction. 
The previous ZEUS measurement [9] was above, but consistent with, the prediction. 

This paper reports a measurement of beauty photoproduction in events with two jets and 
a muon, ep ^ e' bb X ^ e' jjfi X', for < 1 GeV^. 

2 Experimental set-up 

The data sample used in this analysis corresponds to an integrated luminosity C = 110.4± 
2.2 pb"\ collected by the ZEUS detector in the years 1996-1997 and 1999-2000. During 
the 1996-97 data taking, HERA provided collisions between an electron^ beam of Ee = 
27.5 GeV and a proton beam of Ep = 820 GeV, corresponding to a centre-of-mass energy 
^ = 300 GeV (£300 = 38.0 ± 0.6 pb~^). In the years 1999-2000, the proton-beam energy 
was Ep = 920 GeV, corresponding to = 318 GeV (£318 = 72.4 ± 1.6 pb~^). 

A detailed description of the ZEUS detector can be found elsewhere [10]. A brief outline 
of the components that are most relevant for this analysis is given below. 

^ Electrons and positrons arc not distinguished in this paper and are both referred to as electrons. 
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Charged particles are tracked in the central tracking detector (CTD) [11], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organized in nine superlayers covering the 
polar-angle^ region 15° < ^ < 164°. The transverse-momentum resolution for full-length 
tracks is (t{pt)/pt = 0.0058pr © 0.0065 © 0.0014/^^, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [12] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a{E)/E = 0.18/ Ve for electrons 
and a{E)/E = 0.35/-\/E' for hadrons, with E in GeV. 

The muon system consists of rear, barrel (R/BMUON) [13] and forward (FMUON) [10] 
tracking detectors. The B/RMUON consists of limited-streamer (LS) tube chambers 
placed behind the BCAL (RCAL), inside and outside a magnetized iron yoke surrounding 
the CAL. The barrel and rear muon chambers cover polar angles from 34° to 135° and 
from 135° to 171°, respectively. The FMUON consists of six trigger planes of LS tubes 
and four planes of drift chambers covering the angular region from 5° to 32°. The muon 
system exploits the magnetic field of the iron yoke and, in the forward direction, of two 
iron toroids magnetized to ~ 1.6 T to provide an independent measurement of the muon 
momentum. 

The luminosity was measured using the bremsstrahlung process ep epy. The resulting 
small-angle energetic photons were measured by the luminosity monitor [14], a lead- 
scintillator calorimeter placed in the HERA tunnel a.t Z = —107 m. 

3 Data Selection 

The data were selected online by requiring either a high-momentum muon reaching the 
external B/RMUON chambers or two jets reconstructed in the CAL. A dedicated trigger 
requiring two jets and a muon with looser jet and muon thresholds was also used in the 
last part of the data taking. 

Muons were reconstructed offline using the following procedure: a muon track was found 
in the inner and outer B/RMUON chambers or crossing at least 4 FMUON planes, then a 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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match in position and momentum to a CTD track was required. The angular acceptance 
of the F/B/RMUON and of the CTD, and the requirement that the muons reach the 
external chambers define three regions of good acceptance: 

rear -1.6 < r|^' < -0.9, p'^ > 2.5 GeV; 
barrel -0.9 < r]^ < 1.3, pf^ > 2.5 GeV; (1) 
forward 1.48 < r/^ < 2.3, > 4 GeV, pt^ > 1 GeV; 

where ?7^, and pf^ are the muon pseudorapidity, momentum and transverse momentum, 
respectively. 

The hadronic system (including the muon) was reconstructed from energy-flow objects 
(EFOs) [15] which combine the information from calorimetry and tracking, corrected 
for energy loss in dead material. A reconstructed four-momentum [p^xiPyiP^zi "^^s 
assigned to each EFO. 

Jets were reconstructed from EFOs using the kx algorithm [16] in the longitudinally 
invariant mode [17]. The E recombination scheme, which produces massive jets whose 
four-momenta are the sum of the four-momenta of the clustered objects, was used. Muons 
were associated with jets by the kx algorithm: if the EFO corresponding to a reconstructed 
muon was included in a jet then the muon was considered to be associated with the jet. 

The event inelasticity y was reconstructed from the Jacquet-Blondel estimator t/jb = 
{E — pz)/i2Ee) [18], where E — pz = E^ — p^ and the sum runs over all EFOs. 

A sample of events with one muon and two jets was selected by requiring: 

• > 1 muon in one of the three muon-chamber regions defined in Eq. (^; 

• > 2 jets with pseudorapidity |?7'''^*| < 2.5, and transverse momentum p^^ > 7 GeV for 
the highest-p^*^* jet and p^^ > 6 GeV for the second-highest-p^^* jet; 

• that the muon was associated with any jet withp^'^* > 6 GeV and \ri^°^\ < 2.5. To assure 
a reliable p!^^ measurement (see Section the residual jet transverse momentum, 
calculated excluding the associated muon, was required to be greater than 2 GeV; 

• a reconstructed vertex compatible with the nominal interaction point; 

• no scattered-electron candidate found in the GAL; 

• 0.2 < yjB < 0.8. 

The last two cuts suppress background from high-Q^ events and from non-ep interactions, 
and correspond to an effective cut < 1 GeV^ and 0.2 < y < 0.8. 

After this selection, a sample of 3660 events remained. 
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4 Acceptance corrections and background simulation 



To evaluate the detector acceptance and to provide the signal and background distribu- 
tions, Monte Carlo (MC) samples of beauty, charm and light-flavour (LF) events were 
generated, corresponding respectively to six, five and three times the luminosity of the 
data. Pythia 6.2 [19,20] was used as the reference MC, and Herwig 6.1 [21] for 
systematic checks. The branching ratios for direct semi-leptonic b — s> fiX decays and 
for indirect cascade decays into muons via charm, anti-charm, and J/ip, were set to 
-Bdir = 0.106 ± 0.002 and iSindir = 0.103 ± 0.007 [22], respectively. The distribution of 
the decay-lepton momentum in the B-meson centre-of-mass system from Pythia and 
Herwig has been compared with measurements from e~^e~ collisions [23] and found to 
be in good agreement. The generated events were passed through a full simulation of the 
ZEUS detector based on GEANT 3.13 [24]. They were then subjected to the same trigger 
requirements and processed by the same reconstruction programs as for the data. 

Figure H shows the kinematic distributions of p:^^ and 77-'°* for the jet associated with the 
muon, as well as for highest-pT jet that was not associated with a muon (other jet). The 
muon kinematic variables and t]^ are displayed, as well as x}^^^, the fraction of the total 
hadronic E — pz carried by the two highest-pT jets^ 



E-pz 



^jets ^ ^j=l,2y ^ t'Z ) ^ 



The data are compared in shape to the Pythia MC sample in which the relative fractions 
of beauty, charm and LF were mixed according to the cross sections predicted by the 
simulation. The comparison shows that the main features of the dijet-plus-muon sample 
are well reproduced by this MC mixture. The Pythia MC predicts that the non-66 
background comprises 57% prompt muons from charmed-hadron decays and 43% muons 
from light-flavour hadrons, mostly due to in-flight decays of vr and K mesons, with a 
small amount (~ 5% of the LF component) from muons produced in interactions with 
the detector material. The punch-through contribution is negligible. The Herwig Monte 
Carlo (not shown) also gives a good description of the data. 

The detector acceptance for the final cross sections was calculated using the 66 Pythia 
Monte Carlo, in which events were reweighted such that the transverse momentum dis- 
tribution of the 6 quark agreed with that of the NLO calculations. The effects of this 
reweighting on the distributions in Fig^was small. 



^ x^^^'^ is the massive-jets analogue of the x°^^ variable used for massless jets in other ZEUS publications 
[25]. 
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5 Signal extraction and cross section measurement 



Because of the large 6-quark mass, muons from semi-leptonic B-hadron decays tend to 
be produced with large transverse momentum with respect to the direction of the jet 
containing the B-hadron. The beauty signal was extracted by exploiting the distribution 
of the transverse momentum of the muon with respect to the momentum of the rest of 
the associated jet, defined as: 

_ |p"x(pj-^-p^)| 

where p'^ is the muon and pj°* the jet momentum vector. Figure 2a shows the distribu- 
tions, normalized to unity, of the reconstructed muon pJ^' as obtained from the Pythia 
MC, for beauty {ff'^^), charm (/;:^'^^) and LF (/^''■^^) events. The pf distribution 
for beauty peaks at ~ 2 GeV and is well separated from those from charm and LF which 
are peaked close to zero. Since the shapes of /^^''^^ and f^'^^ are very similar, the 
fraction of beauty (a^g) and background (obkg) events in the sample was obtained from a 
two-component fit to the shape of the measured pJ^^ distribution with a beauty and a 
background component: 

/. = Obkg/^'^^ + a,,ff, (4) 

where the p'rf distribution of beauty, ff, was taken from the Pythia MC: ff = fjf'^'^, 
and that of the background, fjf^^, was obtained as explained below. 

The distribution fj^^^ was obtained from the sum of a LF, /^^, and a charm, Z^'^, distri- 
bution weighted according to the charm fraction r obtained from the charm and LF cross 
sections given by Pythia: 

/r'=r/r+(i-o/f- (5) 

The distribution /^^ can be obtained from the distribution of a sample of CTD 
tracks not identified as muons but fulfilling the same momentum and angular requirements 
applied to muons (called "unidentified tracks" in the following). The distribution for 
unidentified tracks, fx, is expected to be similar to /^^, under the assumption that the 
probability for an unidentified track (typically a it ot a. K meson) to be identified as a 
muon, Px^fi, does not depend strongly on pJ^'. This assumption is validated by the MC, 
since the MC distributions for the LF background, /^^'^*", and for the unidentified tracks, 
/^*", are indeed very similar, as shown in Fig. 2a. 

Figure 2b shows the fx distribution obtained from a dijet sample selected without muon 
requirements. The shape obtained from Pythia, /^^, underestimates the tail (for exam- 
ple by 24% at 2.625 GeV). The pJ^' shape of the LF background was therefore obtained 
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as 

xLF.MC 

fLF = f i£ (Q) 

J /J. Jx ; l^; 

where the ratio Z^^'^*"/ is a MC-based correction that accounts for possible differences 
between /^^ and fx due to a residual dependence of Px^^ and to the charm and beauty 
contamination (~ 28% and ~ 2% of respectively) in the dijet sample. 

The data cannot be used to extract the distribution f,';'^. Two cases were therefore con- 
sidered: the distribution given by the Pythia MC, f^'^^-, and the distribution obtained 

rCC.MC 

from the unidentified track sample, as in the case of the LF background: fx jmc ■ The 
average of these two cases was then taken as the nominal f^: 

Figure El shows the result of the pJ^' fit for muons in the rear, barrel and forward regions. 
The sum of the two components reproduces the data reasonably well. The fraction of 
h in the total sample of dijet events with a muon is a^g = 0.224 ± 0.017 (stat.). In the 
determination of the cross sections, the fraction of beauty events in the data was extracted 
by a fit performed in each cross-section bin. 

All the cross sections reported in Section |HJ with the exception of that for h quarks, are 
inclusive muon (or 6-jet) cross sections, obtained by counting muons (or 6-jets) rather 
than events. Muons coming from both direct and indirect h decays are considered to 
be part of the signal. The cross sections are given for dijet events passing the following 
requirements: < 1 GeV^, 0.2 < y < 0.8 and at least two hadron-level jets with 
p^^^ > 7GeV, p^^'^ > 6 GeV and rj , rj ^^^'^ < 2.5. These jets were defined using the 
kx algorithm on stable hadrons, where the weakly decaying B hadrons were considered 
stable. For dijet events with a muon passing the cuts of Eq. ([T)), the acceptance varies 
from 10% at low pj, to 20% at large pj,. 

The cross sections were measured from data collected at two different centre-of-mass 
energies, ^/s = 300 GeV and ^/s = 318 GeV. They were corrected to ^/s = 318 GeV using 
the NLO QCD prediction. The effect of this correction on the final cross section is ~ 2%. 



6 Theoretical predictions and uncertainties 

The measured cross sections are compared to NLO QCD predictions based on the FMNR 
[26] program. The parton distribution functions used for the nominal prediction were 
GRVG-HO [27] for the photon and CTEQ5M [28] for the proton. The 6-quark mass was 
set to rrLf, = 4.75 GeV, and the renormalisation and factorisation scales to the transverse 

mass, fir = f^f= rnT= \ {^P^Y + (Pr)^) + '^ti where p^r^^ is the transverse momentum 

of the h ih) quark in the laboratory frame. Jets were reconstructed by running the /cy 
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algorithm on the four momenta of the b and b quarks and of the third hght parton (if 
present) generated by the program. The fragmentation of the b quark into a B hadron 
was simulated by rescaling the quark three- momentum (in the frame in which = — p^, 
obtained with a boost along Z) according to the Peterson [29] fragmentation function 
with e = 0.0035. The muon momentum was generated isotropically in the B-hadron 
rest frame from the decay spectrum given by Pythia, which is in good agreement with 
measurements made at B factories [23]. 

To evaluate the uncertainty on the NLO calculations, the fe-quark mass and the renor- 
malisation and factorisation scales were varied simultaneously, to maximize the change, 
from rUb = 4.5 GeV and fir = fif= 772^/2 to = 5.0 GeV and fir = fif = 2mT, producing 
a variation in the cross section from +34% to —22%. The effect on the cross section 
of a variation of the Peterson parameter e from 0.002 to 0.0055 [30] and of a change of 
the fragmentation function from the Peterson to the Kartvelishvili parametrisation (with 
a = 13, as obtained from comparisons between NLO QCD and MC distributions) [31,32] 
was less than ±3%. The effects of using different sets of parton densities and of a variation 
of the strong coupling constant (Aq^j^ = 0.226 ± 0.025 MeV) were all within ±4%. 

The NLO cross sections, calculated for jets made of partons, were corrected for jet 
hadronization effects to allow a direct comparison with the measured hadron-level cross 
sections. The corrections were derived from the MC simulation as the ratio of the hadron- 
level to the parton-level MC cross section, where the parton level is defined as being the 
result of the parton showering stage of the simulation. The average between the cor- 
rections obtained from Pythia and Herwig was taken as the central value and their 
difference as the uncertainty. The effect of the hadronization correction was largest in the 
rear region, where the cross section was reduced by (20 ±6)% and smallest at large pt^ 
where it was reduced by (3.0 ± 0.3)%. 

The measured cross sections are also compared to the absolute predictions of two MC 
models, Pythia 6.2 and Cascade 1.1. The predictions of Pythia 6.2 were obtained 
[20] by mixing direct- {•yg bb) and resolved-photon {gg bb, qq ^ bb ) flavour- 
creation processes calculated using massive matrix elements and the flavour-excitation 
(FE) processes {bg bg, bq — > bq), in which a heavy quark is extracted from the photon 
or proton parton density. The FE processes contribute about 27% of the total bb cross 
section. The small (~ 5%) contribution from gluon splitting in parton showers {g — >■ bb) 
was not included. The parton density CTEQ4L [33] was used for the proton and GRVG- 
LO [27] for the photon; the 6-quark mass was set to 4.75 GeV and the 6-quark string 
fragmentation was performed according to the Peterson function with e = 0.0041 [34]. 

Cascade [35] is a Monte Carlo implementation of the CCFM evolution equations [36]. 
Heavy-quark production is obtained from the 0{as) matrix elements for the process 'yg* — > 
bb, in which the initial gluon can be off-shell. The gluon density, unintegrated in transverse 
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momentum {kr), is obtained from an analysis of the proton structure functions based on 
the CCFM equations [37]; in the event generation the gluon density used corresponds 
to the set named "J2003 set 2". The mass of the b quark was set to 4.75 GeV and 
was evaluated at the scale m^. As for Pythia, the 6-quark string fragmentation was 
performed according to the Peterson function with e = 0.0041. 

7 Systematic uncertainties 

The main experimental uncertainties are described below: 

• the muon acceptance, including the efficiency of the muon chambers, of the reconstruc- 
tion and of the MUON-CTD matching, is known to about 10% from a study based on 
an independent dimuon sample [38]; 

• the uncertainty on the shape of the LP and charm background was evaluated by: 

— varying the charm fraction in the background, r, by ±20%. This range was ob- 
tained by fixing the absolute charm-MC normalisation to a measurement of the 
charm dijet cross section [39] and using the Pythia or Herwig MC to extrapolate 
to the kinematic range of the present measurement; 

— varying the jet-track association in the unidentified-track sample; 

— extracting /^^ from a sample of unidentified CTD tracks, reweighted with a MC- 
based parametrisation of Px^^ depending on polar angle and momentum; 

— varying the shape of the charm component of the background between the 
prediction from Pythia and the value obtained from the unidentified track sample; 

— using Herwig instead of Pythia to simulate the background. 

The total uncertainty from these sources is about 10%. As a cross-check, a different 
definition of p^' used to extract the beauty fraction, namely the transverse mo- 
mentum of the muon with respect to the whole jet, including the muon itself, as used 
in a previous ZEUS publication [9]. The results were found to be in good agreement; 

• the 2% uncertainty on the direct-decay branching ratio i^dir introduces a 2% uncer- 
tainty on the 6-jet and on the 6-quark cross sections while it has no effect on the visible 
muon cross sections. The 7% uncertainty on the branching ratio for indirect decays 
-Bindir produces an uncertainty of 1% on the measured cross sections; 

• the uncertainties on the dijet selection, on the energy scale, on the jet and yjB resolu- 
tion and trigger efficiency add up to a 7% uncertainty on the cross sections. 

The uncertainty arising from the model dependence of the acceptance was evaluated as 
follows (the effect on the cross sections is shown in parenthesis): 
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• the Peterson fragmentation parameter e in the MC was varied from 0.002 to 0.006 as 
allowed by LEP and SLD measurements [34,40,41]. The Lund-Bowler fragmentation 
function was used as an alternative, both with the default Pythia parameters and 
with parameters taken from OPAL measurements [34] (±2%); 

• instead of using Pythia reweighted to the NLO p\. distribution it was reweighted as 
a function of r]^^^ and to agree with the measured differential distributions (-2%) 
or without reweighting (+2%); 

• the fraction of flavour-excitation events in Pythia was varied up and down by a factor 
two (± 4%), as allowed by comparisons to the x^^^^ distribution of the data; 

• Herwig was used instead of Pythia (-2%). 

The total systematic uncertainty was obtained by adding the above contributions in 
quadrature. A 2% overall normalization uncertainty associated with the luminosity mea- 
surement was not included. 

8 Results 

All the cross sections reported below, except for the 6-quark cross section, are given for 
dijet events with p'f\p'^^^ > 7, 6 GeV, r] i'''^ r] i''^^ < 2.5, < 1 GeV^ and 0.2 < y < 0.8. 

The first set of measurements are beauty cross sections for muons passing the cuts defined 
in Eq. (^J. The results for the forward, barrel and rear muon-chamber regions are shown 
in Fig. 13^ and Tableland compared with the NLO prediction and MC models. Both the 
NLO and the MC models are in reasonable agreement with the data. 

Figure Eb and Table El show the differential cross-section da/dx}^^^ for muons in the range 
defined by Eq. ((H) in dijet events. The x^°^^ variable corresponds at leading order to the 
fraction of the exchanged-photon momentum in the hard scattering process. It provides 
a tool to measure the relative importance of photon-gluon fusion, 'yg — > bb, which gives 
a peak at x^^^^ ~ 1, and of other contributions, such as gluon-gluon fusion (with a gluon 
coming from the photon) or higher-order diagrams, which are distributed over the whole 
^jets j-ange. The sample is dominated by the high-xJ^''*'^ peak but a low-x^^^*^ component 
is also apparent. The NLO QCD prediction describes the distribution well. Pythia 
is also able to give a good description of the data due to the large contribution from 
flavour excitation at low x^^^'^. Cascade, which generates low-xi^®*'' events via initial-state 
radiation without using a parton density in the photon, tends to underestimate the cross 
section at low x^^^^. 

The differential cross sections in the muonic variables were measured for pt^ > 2.5 GeV 
and —1.6 < r]^ < 2.3. Figure El and Table El show the differential cross-sections da/drj^ 
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and da/dpi^ for muons in dijet events. The NLO QCD predictions and the MC models 
describe the r]^ distribution welL The pt^ distribution is well reproduced by NLO QCD 
while the pi^ slope tends to be too soft in Cascade and Pythia. 

The jet associated with the muon (or /i-jet) reproduces the kinematics of the b (or b) 
quark to a good approximation. The /i-jet is defined as the jet containing the B hadron 
that decays into the muon. Figure ^.-h and Table |3] show the differential cross section 
for the jet associated with a muon passing the cuts of Eq. as a function of the jet 
pseudorapidity, da/dri^~^^^, and transverse momentum, da/dp^^^^, for rj'^'^^^ < 2.5 and 
^M-jet ^ g GeV. The yU-jet distributions are well reproduced by the NLO and by the MC 
models. 

The yU-jet cross sections have been corrected to obtain the cross sections for 6-jets in dijet 
events: a{ep e' jj X). A 6-jet is defined as a jet containing a B (or an anti-B) hadron. 
This correction was performed using Pythia and accounts for the b ^ fi branching ratio, 
including direct and indirect decays, and for the full muon phase space. Figure IHt-d 
and Table El show the differential cross-sections da/drj^^'^^ and da/dp^^'^^. The level of 
agreement of b jets with the NLO QCD and MC predictions is similar to that found for 
the /i jets. It should be noted that the hadronization corrections in the first two rj '^^'^^ 
bins are large (~ —20%). 

To compare the present result with a previous ZEUS measurement given at the fe-quark 
level, the NLO QCD prediction was used to extrapolate the cross section for dijet events 
with a muon to the inclusive 6-quark cross section. The 6-quark differential cross section 
as a function of the quark transverse momentum, da{ep bX)/dp^, for 6-quark pseudo- 
rapidity in the laboratory frame \ri''\ < 2 (corresponding to a rapidity in the ep centre of 
mass of -3.75 < Y^^^ < 0.25), for < 1 GeV^ and 0.2 < y < 0.8, was obtained from the 
dijet cross section for events with a /i-jet within < 2 using the NLO prediction cor- 

rected for hadronization. The b quark was not considered in the definition of the 6-quark 
cross section. As a cross-check, the measurement was corrected to the 6-quark level using 
the Pythia MC, giving a result in agreement within 6%. The result, shown in Figure 
|71 and Table El is compared to the previous ZEUS measurement [9] of the 6-quark cross 
section for p^ > p™™ = 5 GeV and \r]''\ < 2, translated into a differential cross section 
using the NLO prediction and plotted at the average 6-quark transverse momentum, (p^), 
of the accepted events taken from the Monte Carlo: 



The two independent measurements are consistent and in agreement with the NLO QCD 
predictions. 
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9 Conclusions 



Beauty production identified by semi-leptonic decay into muons has been measured in 
dijet events with < 1 GeV^. Differential cross sections for the reaction ep ^ e' bb X 
e' jjfi X' have been measured as a function of the pseudorapidity and transverse momen- 
tum of the muon and of x^°^^. Differential cross sections for the production of 6-jets were 
also measured. 

The results were compared to MC models and to a NLO QCD prediction combined with 
fragmentation and B-hadron decay models. This prediction is in good agreement with 
the data in all the differential distributions. The Pythia MC model is also able to give 
a reasonable description of the differential cross sections. The Cascade MC model also 
gives a reasonable description of the data, except for the low-xi^®*^ region. 

The large excess of the first measurement of beauty photoproduction over NLO QCD, 
reported by the HI collaboration [8], is not confirmed. The present result is consistent 
with the previous ZEUS measurement using semi-leptonic B decays into electrons [9]. 
Beauty photoproduction in ep collisions is reasonably well described both by NLO QCD 
and by a MC model that includes a substantial flavour excitation component. 
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/i-chambers 


Muons 




a ± stat. ± syst. 
(pb) 


(pb) 


Chad 


rear 


484 


0.15 


6.5± 1.5^i;° 


^•'-'-1.0 


0.80 


barrel 


2316 


0.25 


38.2 ± 3Atli 


oo q+11-0 


0.89 


forward 


868 


0.21 


16.6±3.3j:|^ 


6.51?:^ 


0.86 



Table 1: For each muon-chamher region defined in Eq. (Qp the columns show: 
the number of selected muons; the beauty fraction a^^ obtained from the p^j?' fit; the 
measured beauty cross section with the statistical and systematic uncertainties; the 
NLO QCD prediction corrected to the hadron level with the theoretical uncertainty 
and the hadronization correction. For further details see the caption to Fig. ^ 



^jetS j^g^j^gg 


da/dx^^^^ ± stat. ± syst. 
(pb) 


Chad 


0.00, 0.25 


12.9 ± 6.7^^:^ 


0.68 


0.25, 0.50 


20.8 ± 7.2+?:^ 


0.83 


0.50, 0.75 


30.0 ± 5.2+^:^ 


0.86 


0.75, 1.00 


165 ± 14+22 


0.90 



Table 2: Differential muon cross section as a function of x^^^^ . The multiplicative 
hadronization correction applied to the NLO prediciton is shown in the last column. 
For further details see the caption to Fig. ^ 



rj^ range 


da /drj^ ± stat. ± syst. 
(pb) 


Chad 


-1.6, -0.75 


4.6±1.3+°| 


0.83 


-0.75, 0.25 


18.8±2.2l|^ 


0.88 


0.25, 1.30 


16.7±2.4l|^ 


0.92 


1.30, 2.30 


10.0±2.3l^:^ 


0.91 


range 
(GeV) 


da/dpf^ ± stat. ± syst. 
(pb/GeV) 


Chad 


2.5, 4.0 


16.1 ±2.3+1? 


0.87 


4.0, 6.0 


7.1 ± i.otli 


0.92 


6.0, 10.0 


1.69±0.32+°:22 


0.97 



Table 3: Differential muon cross section as a function oft]^ andpi^. For further 
details see the caption to Fig. 
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^fi jet raj^gg 


da/drj'^ j''* ± stat. ± syst. 
(pb) 


Chad 


-1.6,-0.6 


8.6±1.7t\i 


0.73 


-0.6, 0.4 


18.9 ± 2.2+^-8 


0.84 


0.4, 1.4 


14.3 ± 2.3tli 


0.97 


1.4, 2.5 


UA ± 2.7tli 


1.00 


plj, ■' range 


da/dpif^^^ ± stat. ± syst. 


Chad 


(GeV) 


(pb/GeV) 




6,11 


6.41 ± 0.67+^;°° 


0.88 


11,16 


2.98 ± 0.37+°;| 


0.88 


16,30 


0.51 ±0.1 iiJ]:J]^ 


0.90 



Table 4: Differential cross section for jets associated with a muon as a function 
qJ ^f^-ict andp^^^'^^. For further details see the caption to Fig.\^ 





da/d-q^^""^ ± stat. ± syst. 
(pb) 


Chad 


-1.6,-0.6 


152 ± 29tl\ 


0.68 


-0.6, 0.4 


356 ± 4ll| 


0.78 


0.4, 1.4 


275 ± 45l| 


0.96 


1.4, 2.5 


229 ± AAtfsi 


1.06 


range 


da /dp^^^^ ± stat. ± syst. 


Chad 


(GeV) 


(pb/GeV) 




6,11 


137 ± 141^1 


0.85 


11,16 


43.8 ± 5.51^2^0 


0.86 


16,30 


5.7±1.2;j;° 


0.89 



Table 5: Differential cross section for h-jets as a function of rf j'^* andp\.^°^ . For 
further details see the caption to Fig. O 
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u-ict 

p>P range 


da/dpt^-^^^ ± stat. ± syst. 


Chad 




da/dp^ ± stat. ± syst. 


(GeV) 


(|^M-jct| <- 2) (pb/GeV) 




(GeV) 


(pb/GeV) 


6,11 


6.22 ± O.QStiM 


0.87 


8.5 


90 ± 9^}^ 


11,16 


2.74 ± OMtof^ 


0.88 


13.6 


18.2 ± 2.41^:^ 


16,30 


o.43±o.ioiEJ:[!^ 


0.88 


21.25 


2.0±0.5l^j 



Table 6: Differential b-quark cross-section da/dp^. The first columns show the 
differential fi-jet cross-section da/dpt^-^'^^ restricted to < 2, the corresponding 

hadronization correction and the average b-quark transverse momentum, (p^), as 
obtained from the Pythia MC. The b-quark differential cross section da/dp^ for 
|?7^| < 2, evaluated at {p^), is shown in the last column. For further details see the 
caption to Fig. 
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Figure 1: Distributions for the dijet-plus-muon sample (points) compared to the 
predictions of the Pythia Monte Carlo (full line) normalized to the data. The 
shaded histogram shows the beauty component and the dotted line is the sum of 
charm and beauty. The plots show (a) the transverse momentum of the jet asso- 
ciated with the muon; (b) its pseudorapidity; (c) the transverse momentum of the 
highest-pT non-^- associated jet; (d) its pseudorapidity; (e) the transverse momen- 
tum of the muon; (f) its pseudorapidity; (g) the distribution of x^^^^ . 
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Figure 2: (a) The pJ^^ distribution as predicted by the Pythia Monte Carlo for 
reconstructed muons from beauty (f^^'^^ , dashed histogram), charm (f'^'^'^, dotted 
histogram) and light- flavours (fjf'^^, dash-dotted histogram), and for unidentified 
tracks (f^^ , full-line histogram). The distributions are normalized to unity, (b) 
Thep^rf distribution of unidentified data tracks (points), compared to the prediction 
from Pythia (full line). The charm and beauty components are also shown as the 
dotted- and dashed-line histograms, respectively. 
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Figure 3: The distribution for events with muons in the rear, barrel and 
forward regions defined in Eq. (OP- The data (points) are compared to the mixture 
of beauty and charm+LF background obtained from the p!^^ fit (full line). The 
shaded histogram represents the beauty component while the dashed-line histogram 
is the background. 
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Figure 4: Cross sections for muons coming from b decays in dijet events with 
p'^^' > 7,p^T*' > 6 GeV, rf^^^rf^^^ < 2.5, 0.2 < ?/ < 0.8, < i Q^y^ passing the 
selection of Eq. 0). (a) The cross section for the forward, barrel and rear regions 
(defined in Eq. (b) The differential cross section as a function of x^^^^ . The 

data (points) are compared to the predictions of NLO QCD (dotted line: parton- 
leveljets; dashed line: jets corrected to the hadron level). The full error bars are the 
quadratic sum of the statistical (inner part) and systematic uncertainties. The band 
around the NLO prediction represents the variation on the theoretical predictions 
obtained by varying the b-quark mass, fir and fif, as explained in the text. The 
data are also compared to the predictions of the Pythia (solid line histogram) and 
Cascade (dot-dashed line histogram) Monte Carlo models. 
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Figure 5: Differential cross sections as a function of (a) the muon pseudorapidity 
Tj^ and (b) transverse momentum p^, for pf^ > 2.5 GeV and —1.6 < rj^ < 2.3, 
for muons coming from b decays in dijet events with p^^^^ > 7,p!^^'^ > QGeV, 
j^jcti ^ j^jct2 ^ 2.5, 0.2 < y < 0.8, < 1 GeV"^ . The data (points) are compared to 
the predictions of NLO QGD (dotted line: parton-level jets; dashed line: corrected 
to hadron-level jets). The full error bars are the quadratic sum of the statistical 
(inner part) and systematic uncertainties. The band around the NLO prediction 
represents the variation on the theoretical predictions obtained by varying the b- 
quark mass and fir and fif as explained in the text. The data are also compared to 
the predictions of the Pythia (solid line histogram) and CASCADE (dot- dashed line 
histogram) Monte Garlo models. 
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Figure 6: Differential cross sections as a function of (a) the pseudorapidity 
j^M-jet g^^^ transverse momentum pf^~^'^^ of the jet associated to the muon, for 

> 6 GeV, r/'^-j'=* < 2.5, for muons passing the selection of Eq. ^ and coming 
from b decays; differential cross sections as a function of (c) the pseudorapidity 
^b-jet transverse momentum p^"'''^* of the jet containing a B hadron. 

All the cross sections are evaluated for dijet events with p^^*^ > 7,pJ^*^ > 6 GeV , 
^jeti^j^jeta ^ 2.5, 0.2 < y < 0.8, Q'^ < 1 GeV"^. The data (points) are compared to 
the predictions of NLO QCD (dotted line: parton level; dashed line: corrected to 
hadron level). The full error bars are the quadratic sum of the statistical (inner part) 
and systematic uncertainties. The band around the NLO prediction represents the 
uncertainty on the theoretical prediction corrected for hadronization. The data are 
also compared to the predictions of the Pythia (solid line histogram) and Cascade 
(dot-dashed line histogram) Monte Carlo models. 
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Figure 7: Differential cross section for b-quark production as a function of the 
b-quark transverse momentum for b-quark pseudorapidity \ri''\ < 2 and for < 
1 GeV'^ , 0.2 < y < 0.8. The filled points show the ZEUS results from this analysis 
and the open point is the previous ZEUS measurement in the electron channel 
[9]. The full error bars are the quadratic sum of the statistical (inner part) and 
systematic uncertainties. The dashed line shows the NLO QCD prediction with the 
theoretical uncertainty shown as the shaded band. 
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